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Abstract 



to vector couplings -r— 
JV 



The mass mjf« and vector coupling fx* of the i^*-meson, as well as the ratio of the tensor 
fx 

, are computed in lattice QCD. Our simulations are performed 

K* 

in a partially quenched setup, with two dynamical (sea) Wilson quark flavours, having a 
maximally twisted mass term. Valence quarks are either of the standard or the Osterwalder- 
Seiler maximally twisted variety. Results obtained at three values of the lattice spacing are 
extrapolated to the continuum, giving mx* = 981(33)MeV, fx* = 240(18)MeV and 

^ 2GeV ) - 0.704(41). 



fv 



1 Basics 



The aim of the present letter is to present novel lattice results for the mass of the K*- 
meson, as well as its vector and tensor couplings (/y and fx respectively), defined 
in Euclidean space-time as follows: 



(0\Vj\K*;\) = -if V e$m K * , (1.1) 
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(0\T Oj \K*;X) = -if T e]m K * . (1.2) 

In the above expressions, Vj = sjjd is the vector current (spatial components only; 
j = 1,2,3), Tqj = isaojd is the tensor bilinear operator (temporal component), and 
6j denotes the polarization vector. 

Our results are based on simulations of the ETM Collaboration (ETMC) pQ, 
with Nj = 2 dynamical flavours (sea quarks) and "lightish" pseudoscalar meson 
masses in the range 280 MeV < mps < 550 MeV. With three lattice spacings (a = 
0.065 fm, 0.085 fm and 0.1 fm) we are able to extrapolate our results to the continuum 
limit. Our simulations are performed with the tree-level Symanzik improved gauge 
action. For the quark fields we adopt a somewhat different regularization for sea 
and valence quarks. The sea quark lattice action is the so-called maximally twisted 
standard tmQCD (referred to as "standard tmQCD case") [2]. The Nj = 2 light sea 
quark flavours form a flavour doublet x = (« , d) and the fermion lattice Lagrangian 
in the so-called "twisted basis" is given by 



Am = X 



D w + ifig^r 3 X , (1-3) 



where r 3 is the isospin Pauli matrix and D\y denotes the critical Wilson-Dirac 
operator. By "critical" we mean that, besides the standard kinetic and Wilson 
terms, the operator also includes a standard, non-twisted mass term, tuned at the 
critical value of the quark mass (k ct in the language of the hopping parameter), so 
as to ensure maximal twist. With only two light dynamical flavours, strangeness 
clearly enters the game in a partially quenched context. For the valence quarks we 
use the so-called Osterwalder-Seiler variant of tmQCD, which consists in maximally 
twisted flavours which, unlike the standard tmQCD CclSG ; £11*6 not combined into 
isospin doublets: 

£os = ^2 Qf\Dw +*A*/75 Qf ! (1-4) 

f=d,s 

with sign(/if) = ±1 (see below for details). This action, introduced in ref. [3] and 
implicitly used in [3], has been studied in detail in ref. [5]. For the case in hand 
(i.e. if *-related quantities) we only need down- and strange-quark flavours in the 
valence sector. Note that the choice of maximally twisted sea and valence quarks 
implies 0(a)-improvement of the physical quantities (i.e. the so-called automatic 
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improvement of masses, correlation functions and matrix elements) [6]. Thus unitar- 
ity violation, which plagues any partially quenched theory at finite lattice spacing, 
is an 0(a 2 ) effect. 

The sign of fj, s may be that of \x& or its opposite. We conventionally refer to the 
setup in which sign(^d) = — sign(// s ) as the "standard twisted mass regularization" 
(denoted by tm) and the setup with sign(^d) = sign(/i s ) as the "Osterwalder-Seiler 
regularization" (denoted as OS). Quenched pseudoscalar masses and decay constants 
in tm- and OS-setups have already been studied [HE]. 

The continuum operators of interest are expressed, in terms of their lattice 
counterparts, as follows: 



t Acont 
/y-icont 



ZaA? 



+ 0{a 2 ) 
+ 0(a 2 ) 



= z v v™ 



+ 0(a 2 



0(< 



a 



1.5) 
1.6) 



where T, 



The vector and axial currents are normalized by the scale 
independent factors Zy and Za, while Zt = Zt{^) runs with a renormalization 
scale fj, (i.e. it is defined in a given renormalization scheme). 

The vector boson mass, my, as well as fy and /r, are obtained form two-point 
correlation functions at zero spatial momenta and large time separations. These are 
defined in the continuum (Euclidean space-time) as 



\x° 



f V my 



-y. 

3 Zs 

j 

2 



exp[— myT/2] cosh [my(T/2 
d 3 x (T 0j (x) 4(0)) c 



*°)] 



\ cont 



(1.7) 



exp[-m v T/2] cosh [my (T/2 - x )] . (li 



The asymptotic expressions of the above equations correspond to the large time 
limit of the correlation functions (symmetrized in time), with periodic boundary 
conditions for the gauge fields and (anti)periodic ones for the fermion fields in the 
(time)space directions (i.e. <C x° <C T/2). These are actually the boundary 
conditions of our lattice simulations. The lattice correlation functions are related to 
the continuum ones as suggested by eqs. (|1.5p . (|1.6|) : 



Cy" nt (x°) 

Cp nt (x°) 



Zq-< C' 



tm/ 



T 



X 



+ 0(a 2 ) 
+ 0{a 2 ) 



Z 2 yC^° 

z\c% 



OS 

f 



+ 0{a 2 ) 
(x ) + 0{a 2 ) 



(1.9) 
(1.10) 



The meaning of the notation C A m , Cj~ , etc. should be transparent to the reader. 
The ratio Jt / fv is computed from the square root of the ratio of correlations func- 
tions C^j? nt /Cy nt , in which many systematic effects cancel. We compute the vector 
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meson mass and decay constant from Cy nt and the ratio fr/fv from the ratio 
of correlation functions C£ ont /C£? nt . The tensor coupling fx is then obtained by 
multiplying fx/fv by fv- 

Note that fv is a scale independent quantity, while fx (/J-) depends on the renor- 
malization scale fj,, as well as the renormalization scheme. The scale and scheme 
dependence of the latter quantity is carried by the renormalization factor Zt([a); we 
opt for the MS-scheme and for fi = 2 GeV. 

2 Results 

ETMC has generated Nf = 2 configuration ensembles at four values of the inverse 
gauge coupling; in this work we make use of only three of them. Light mesons consist 
of a valence quark doublet, with twisted mass a/i? equal to that of the sea quarks; 
afi£ = a/x sea . Heavy-light mesons consist of a valence quark pair {apu = a/x sea , afih)- 
As already stated, these bare quark mass parameters are chosen so as to have light 
pseudoscalar mesons ("pions") in the range of 280 < mpg < 550 MeV and heavy- 



light pseudoscalar 


mesons ("Kaons' 


') in the 


range 450 


< rnps 


< 650 MeV. The 


simulation parameters are gathered 


in Table [TJ 








P 


a" 4 (L 3 x T) 


afi£ = 




afih 




^meas 


3.80 


2\ 6 x 48 


0.0080, 


0.0110 


0.0165, 


0.0200 


180 


(o ~ 0.1 fm) 








0.0250 






3.90 


24 3 x 48 


0.0040 




0.0150, 


0.0220 


400 










0.0270 








24 3 x 48 


0.0064, 


0.0085, 


0.0150, 


0.0220 


200 






0.0100 




0.0270 






3.90 


32 3 x 64 


0.0030, 


0.0040 


0.0150, 


0.0220 


270/170 


(a ~ 0.085 fm) 








0.0270 






4.05 


32 3 x 64 


0.0030, 


0.0060, 


0.0120, 


0.0150 


200 


(a ~ 0.065 fm) 




0.0080 




0.0180 







Table 1: Simulation details 

Our calibrations are based on earlier collaboration results. The ratio ro/o, 
known at each value of the gauge coupling f3 from ref. [9] , allows to express our raw 
dimensionless data (quark masses, meson masses and decay constants) in units of 
tq. Knowledge of the renormalization constant Zp in the MS scheme at 2 GeV (see 
ref. |10j ) enables us to pass from bare quark masses to renormalized ones (again in ro 
units). Using only data with light valence quarks in the tm-setup, we have applied 
the procedure described in refs. [T][9] for the determination of the physical continuum 
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light quark mass fJ>^h- From the data concerning light and heavy valence quark 
masses in the tm-setup f9] , we determine the physical continuum strange quark mass 
//^ s (2GeV). These quark mass values are listed in Table O The Sommer scale we 
use, based on an analysis with three values of the lattice spacing, is ro = 0.448(5) fm. 
This updates our previous ro computation, derived with two /3's, cf. ref. [I]. 

We see from eqs. (|1.9|) and (|1,1U|) that we need to know the renormalization 
parameters Zy, Za, and Zt- These quantities, as well as Zp, have been computed 
in ref. [ID] , in the RI/MOM scheme; Zp and Zt are perturbatively converted to MS. 
In the same work a Zy estimate, obtained from a Ward identity, is also provided. 
In Table [2] we gather the most reliable estimates of ref. [TU], which we have used in 
the present analysis, as well as our estimates of the r^/a ratio. 



p 


Zy 


Z A 


Z T ls {2 GeV) 


Z P ls (2 GeV) 


r /a 


3.80 


0.5816(02) 


0.746(11) 


0.733(09) 


0.411(12) 


4.54(07) 


3.90 


0.6103(03) 


0.746(06) 


0.743(05) 


0.437(07) 


5.35(04) 


4.05 


0.6451(03) 


0.772(06) 


0.777(06) 


0.477(06) 


6.71(04) 



Table 2: The renormalization parameters used in our analysis and the ro/a values at 
each gauge coupling. Zy is obtained from a lattice vector Ward identity, while the 
other renormalization constants are obtained from the RI/MOM scheme; for details 
see ref. [TO] . 



^fp GeV) ^ s (2 GeV) 



3.6(2) MeV 95(6) MeV 



Table 3: The quark mass values (in the MS scheme), used in our analysis; see ref. [9]. 

As can be seen in Table [H at (3 = 3.90 we have performed more extensive 
simulations, which enable us to check in some detail the quality and stability of the 
measured physical quantities. We wish to highlight straightaway the two problems 
we have encountered in these tests, performed for the tm-setup: (i) For all sea quark 
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masses, when the valence quark attains its lightest value a[n = 0.0040, the vector 
meson effective mass has a poor plateau. The situation already improves at the 
next quark mass a/j,£ = 0.0064. Nevertheless, since the signal-to-noise ratio behaves 
as expected (i.e. it drops like exp[— (my — mps)x ]) the p-meson mass and decay 
constant can still be extracted (see results presented in ref. [H]). (h) A poor quality 
vector meson effective mass is also seen when < ^ S ea- This problem is absent in 
the pseudoscalar channel. 

The above problems are easily avoided in the present work, since the quantities 
of interest are related to the i^T*-meson, consisting of a down and a strange valence 
quark mass (/J- u /d < Us)- We thus proceed as follows: at each f3 value, we compute 
the necessary observables (vector meson mass my, vector decay constant fy, and 
the ratio fa/fv), for all combinations of a\it = a^ sea and afih (with m < fi^). In 
this way unitarity holds in the light quark sector, while the heavy valence quark 
mass, in a partially quenched rationale, spans a range around the physical value /i s . 
Examples of the quality of our signal are given in Figs. [1] and [2j the lightest mass 
is a/imin and the heavy mass, corresponding to the physical strange value a^ s , is 
obtained by interpolation, as will be explained below. 
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Figure 1: Effective vector meson mass r^my at three values of the lattice spacing. 
The light quark mass is o/i m i n (see Table [Q) and the heavy quark mass a/x^ is close 
to that of the physical strange quark, (a) tm-setup ; (b) OS-setup. Plateau intervals 
are indicated by straight lines. 



Statistical errors are estimated with the bootstrap method, employing 1000 
bootstrap samples. A reliable direct determination of the ratio fa/ fv i n the OS- 
setup is not possible, because the ratio of correlation functions C? s /Cy S do not 
display satisfactory plateaux, due to big statistical fluctuations of the tensor cor- 
relator C? s . We only present fa/fv results in the tm-setup, obtained from the 
better-behaved correlation function Cf 11 . In Fig. [3] we show results for this ratio at 
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a/i m ; n and also at a heavier light quark mass. 
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Figure 2: Vector decay constant r^fy at three values of the lattice spacing. The 
light quark mass is a// m i n (see Table [1]) and the heavy quark mass a\Xh is close to 
that of the physical strange quark, (a) tm-setup ; (b) OS-setup. Plateaux intervals 
are indicated by straight lines. 
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Figure 3: The ratio fr/fv i n the tm-setup, at three values of the lattice spacing 
and heavy quark mass afih, close to that of the physical strange quark, (a) For 
the lightest quark mass a/i m i n ; (b) for the next-to-lightest quark mass. Plateaux 
intervals are indicated by straight lines. 



Regarding vector meson masses my and couplings fy, both tm- and OS-results 
display similar plateau quality and statistical accuracy. At finite lattice spacing and 
for equal bare quark masses, tm- and OS-estimates of my are compatible within 
errors. Agreement is also very good for fy, with occasional discrepancies, interpreted 
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as cutoff effects, showing up at the coarsest latticqj. Contrary to the well known 
large 0(a?) isospin breaking effects in the neutral to charged pion splitting mass, 
no numerically large differences are observed between tm and OS results for fy and 
my. This fact is in agreement with theoretical expectations, see ref. |12j . 



P m r m^(ls) r m» s (f, S ) r f^(£,s) r f° s (e,s) [f T /fv] tm (ls) 



3.80 


0.0080 




2.443(41) 


2.471(30) 


0.642(18) 


0.700(13) 


0.764(38) 




0.0110 




2.508(32) 


2.500(23) 


0.651(14) 


0.706(15) 


0.792(35) 


3.90 


0.0040 




2.410(41) 


2.381(38) 


0.610(21) 


0.643(17) 


0.755(19) 




0.0064 




2.441(32) 


2.427(35) 


0.626(22) 


0.659(12) 


0.726(20) 




0.0085 




2.484(48) 


2.441(33) 


0.628(16) 


0.652(16) 


0.776(27) 




0.0100 




2.468(54) 


2.481(32) 


0.619(20) 


0.657(16) 


0.774(31) 




0.0030(L= 


=32) 


2.259(75) 


2.335(45) 


0.577(20) 


0.639(16) 


0.714(20) 




0.0040(L= 


=32) 


2.364(32) 


2.371(50) 


0.599(22) 


0.640(21) 


0.722(19) 


4.05 


0.0030 




2.305(86) 


2.263(80) 


0.568(49) 


0.588(40) 


0.742(27) 




0.0060 




2.439(67) 


2.295(76) 


0.618(41) 


0.578(46) 


0.768(30) 




0.0080 




2.512(65) 


2.427(48) 


0.649(31) 


0.648(27) 


0.741(31) 


CL 


t- L u/d 




2.227(71) 


2.200(60) 


0.545(41) 


0.525(30) 


0.701(46) 


expt. 






2.025 




0.493 







Table 4: Results for three values of lattice spacing and several light quark masses 
afi£, interpolated to the physical strange mass fj, s . Vector mass and vector decay 
constant results are presented for both tm- and OS-setups . The ratio /t/Jv results 
are given only in the tm-setup. Our extrapolations at the \i u u physical point and 
in the continuum limit are also shown. In the last row the experimental results for 
the vector mass and the vector decay constant, in units of ro, have been added. 

The extrapolation to the physical quark masses is carried out in two steps. 
First, for fixed values of the gauge coupling (3 and light quark mass afj,£ = a^ sea , 
we perform linear interpolations of romy, r§fy and /r//y to the physical strange 
quark mass /j, s . The second step consists in using these interpolated results for a 
combined fit of our data at three lattice spacings and all available light quark masses, 
in order to determine the continuum value of the quantity of interest (romy, r^fy 
and fr/fv)- The fitting function we use is 

a 2 

m v r Q = C Q (/j, s ro) + Ci(/x s r ) /^r + D(fj, s r Q ) — , (2.1) 

r o 

1 Given the large fluctuations of fr/ fv in the OS-setup at the finer lattice spacing, we only quote 
results for this ratio in the tm-setup. 
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and similarly for fyTo and fr/fv- The results of the interpolations in the heavy 
quark mass /^ to the physical value at each f3 and afi£, are gathered in Tabled! 
In the same Table we also display the results of the combined chiral and continuum 
extrapolations. Note that for the three quantities of interest, mv,fv and fr/Iv, 
the value of x 2 /d.o.f. is less than unity. The linear dependence of our data on the 
light quark mass agrees with the predictions of chiral perturbation theory for the 
ratio fr/ fv hi the K* mass range; see refs. [T3J[T1]. 
Our final results, extracted in the tm-setup, are 

m K . = 981(31)(10)[33]MeV , (2.2) 
f K * = 240(18)(02)[18]MeV . (2.3) 

The first error includes the statistical uncertainty and the systematic effects re- 
lated to the simultaneous chiral and continuum fits, mass interpolations and ex- 
trapolations, and uncertainties in the renormalization parameters. The second error 
arises from that of tq. These two errors, combined in quadrature, give the to- 
tal error in the square brackets. It is encouraging that these results agree with 
the ones obtained in the OS-setup (which is a different regularization) , namely 
m K * = 969(27) (10) [29]MeV and f K * = 231(13)(02)[13]MeV. Compared to the ex- 
perimentally known values, txik* = 892MeV and fK* = 217MeV, the vector meson 
mass is 2-3 standard deviations off, while the decay constant is compatible within 
about one standard deviation. 

Our final estimate (tm-setup) for the ratio of vector meson couplings is 



/ T (2GeV) 



fv 



= 0.704(41) . (2.4) 

K* 



This is compatible with the continuum limit quenched result [f T (2GeV)/ fv]K* = 
0.739(17)(3) of ref. [15]. We are also in agreement with the result of the RBC/UKQCD 
collaboration [IB]; using Nf = 2 + 1 dynamical fermions at a single lattice spacing, 
they quote [/t(2 GeV)//v]jf* = 0.712(12). The lattice results are also in agreement 
with the sum rules' estimate [/t(2 GeV)/fv]K* = 0.73(4), quoted in [17] , 
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Figure 4: r^my plotted against the renormalized light quark mass ro/tf ; (a) tm- 
setup; (b) OS-setup. The continuous lines are combined chiral and continuum ex- 
trapolations to the physical point. The bottom (black) line corresponds to eq. (|2,1|) 
at a = 0. The separation among the four lines in (a) is invisible to the naked eye 
(i.e. small scaling violations). 




Figure 5: rofv plotted against the renormalized light quark mass rofi£ ; (a) tm-setup; 
(b) OS-setup. The continuous lines are combined chiral and continuum extrapola- 
tions to the physical point. The bottom (black) line corresponds to eq. (|2.ip at 
a = 0. 
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